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Neuroimaging studies have demonstrated that the phenomenology of schizophrenia maps onto diffuse alter- 
ations in large-scale functional and structural brain networks. However, the relationship between structural 
and functional deficits remains unclear. To answer this question, patients with established schizophrenia and 
matched healthy controls underwent resting-state functional and diffusion weighted imaging. The network- 
based statistic was used to characterize between-group differences in whole-brain functional connectivity. 
Indices of white matter integrity were then estimated to assess the structural correlates of the functional al- 
terations observed in patients. Finally, group differences in the relationship between indices of functional and 
structural brain connectivity were determined. Compared to controls, patients with schizophrenia showed 
decreased functional connectivity and impaired white matter integrity in a distributed network encom- 
passing frontal, temporal, thalamic, and striatal regions. In controls, strong interregional coupling in neural 
activity was associated with well-myelinated white matter pathways in this network. This correspondence 
between structure and function appeared to be absent in patients with schizophrenia. In two additional 
disrupted functional networks, encompassing parietal, occipital, and temporal cortices, the relationship be- 
tween function and structure was not affected. Overall, results from this study highlight the importance of 
considering not only the separable impact of functional and structural connectivity deficits on the pathoaeti- 
ology of schizophrenia, but also the implications of the complex nature of their interaction. More specifically, 
our findings support the core nature of fronto-striatal, fronto-thalamic, and fronto-temporal abnormalities 
in the schizophrenia connectome. 

© 2014 The Authors. Published by Elsevier Inc. 
This is an open access article under the CC BY-NC-ND license 
(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Schizophrenia is a debilitating mental illness characterized by 
complex cognitive, affective, and behavioural symptoms that signif- 
icantly deteriorate an individual's quality of life (Hack and Newhill, 
2007). Studies adopting magnetic resonance imaging have shown that 
the complexity and severity of these various symptoms are mirrored 
in the extent of focal grey and white matter alterations in the brain 
(Bora et al., 201 1 ; Fornito et al., 2009; Shepherd et al., 2012). A recent 
review (Fornito et al., 2012) has also drawn attention to the repro- 
ducible and widespread deficits in functional (Alexander-Bloch et al., 
2010; Bassett et al., 2012; Bassett et al., 2008; Fornito et al., 2011; 
Zalesky et al., 201 1 ; Lynall et al., 2010; Rubinov et al., 2009; Liu et al., 
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2008; Collin and Pol, 2011; Wang et al., 2010) and structural brain 
connectivity (Zalesky et al., 2011; Collin et al., 2013; Pettersson-Yeo 
et al., 2011) that are characteristic of schizophrenia. These connec- 
tivity abnormalities support the conceptualization of schizophrenia 
as a disconnection syndrome (Friston, 2005; Friston and Frith, 1995). 
Accordingly, a comprehensive characterization of the schizophrenia 
connectome is an important step to understanding the neural basis of 
this disorder and to developing more accurate biomarkers for diag- 
nosis and treatment efficacy (Bassett et al., 2008; Cocchi et al., 2012; 
Insel and Wang, 2010; Fornito and Harrison, 2012; van den Heuvel et 
al.,2010; Filippi et al., 2013). 

The functioning of the human brain is constrained by its anatomi- 
cal interconnections (Bassett et al., 2008; van den Heuvel et al., 2010; 
Sporns, 2011; Sporns, 2012; Honey et al., 2009; Hermundstad et al., 
2013; Deco et al., 2013; Skudlarski et al., 2008). The most striking 
example of this can be found in the healthy brain, where the pres- 
ence of well-myelinated white matter pathways is known to be an 
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accurate predictor of synchronization in neural activity measured be- 
tween corresponding regions (i.e., functional connectivity, Honey et 
al., 2009). The converse is not necessarily true; that is, functional 
connectivity can be established in the absence of a direct structural 
pathway (Stephan et al., 2009; Stephan et al., 2006). For example, the 
functional interplay between two brain regions may be mediated by 
a third region. Whether the structure-function connectivity relation- 
ship is preserved in schizophrenia and in other psychiatric illnesses 
remains a largely unresolved matter. 

To date, a small number of neuroimaging studies have reported 
spatial coincidence in aberrant indices of both large-scale functional 
connectivity and structural connectivity in schizophrenia (Yan et al., 
2012; Camchong etal., 2011; Zalesky etal., 2010). For instance, in one 
study concurrently assessing functional and structural connectivity, 
Zalesky et al. (2010) demonstrate overlapping anatomical and func- 
tional alterations in brain networks encompassing the mediofrontal, 
temporal, parietal and occipital regions. Recent studies investigating 
the relation between structural connectivity and functional connec- 
tivity in schizophrenia have reported evidence of both weaker and 
stronger structure-function coupling relative to the healthy brain 
(Skudlarski et al., 2010; van den Heuvel et al., 2013), suggesting that 
variation in this relationship may be cohort and/or network-specific. 
However, it remains unclear whether abnormal structure-function 
relationships are convergent with abnormalities in functional and 
structural connectivity in their own right; such links may point to 
dependencies between metrics of connectomic pathophysiology that 
define this disorder. 

Given the paucity of research in this area, we sought to as- 
sess the structural underpinnings of altered functional connectivity 
in schizophrenia by investigating whether the association between 
functional and structural connectivity was preserved in functional 
networks affected by the disorder. To this end, connectomic tech- 
niques were applied to diffusion and resting-state functional imaging 
data acquired in a sample of well characterized patients with es- 
tablished schizophrenia and a carefully matched group of healthy 
controls. The network-based statistic (NBS, Zalesky et al., 2010) was 
first used to comprehensively test the entire connectome for func- 
tional connections showing evidence of abnormal synchronization 
in patients with schizophrenia. This whole-brain analysis was fol- 
lowed by a selective assessment of the white matter integrity of the 
anatomical pathways interconnecting the pairs of regions identified 
as showing abnormal functional connectivity. Finally, we sought to 
assess whether the relationship between indices of functional and 
structural connectivity in the networks of interest was altered in the 
patient group. 

Based on our previous investigations (Zalesky et al., 2010), we hy- 
pothesize that the schizophrenia connectome comprises disrupted 
white matter pathways that are spatially correspondent with pairs 
of regions showing abnormal resting-state functional connectivity. 
We further hypothesize that altered functional and structural net- 
work connectivity may be linked to a loss of the normal coherence 
between these measures observed in the healthy brain. Specifically, 
we expect that stronger structure-function correlations in patients 
with schizophrenia compared to controls may be a marker of greater 
rigidity in functional network dynamics of schizophrenia. Conversely, 
weaker coupling may indicate a randomization of the normal con- 
straints anatomical interconnections impart on brain function. 

2. Materials and methods 

2.3. Participants 

Eighteen adults meeting the diagnostic criteria for schizophrenia 
(Structured Clinical Interview for DSM-lVAxis-1 Disorders, First, 1997) 
and 18 healthy controls were recruited from the general community. 
Controls did not meet lifetime criteria for any Axis-1 disorder or have 



a family history of psychosis. Patients with comorbid Axis-I disor- 
ders were also excluded, excepting past major depressive episodes 
or past substance abuse/dependence. Additional exclusion criteria 
included neurological conditions, significant head injury, electrocon- 
vulsive therapy, systemic steroid use, intellectual disability (IQ. < 80), 
MRl incompatibility or intolerance (e.g. cardiac pace-maker), acute 
intoxication, and excessive head motion (>2 mm) during fMRI data 
acquisition. 

One patient was excluded due to excessive head motion and one 
control participant due to scanner intolerance. The remaining 1 7 par- 
ticipants in each cohort were closely matched for age, gender, years of 
education, handedness, and both current and historical substance use 
(with the exception of greater tobacco intake in patients. Table 1 ). As 
is generally reported, the schizophrenia cohort had a lower full-scale 
lQ(Wechsler, 1999). 

All patients with schizophrenia were compliant to an atypical 
neuroleptic pharmacotherapy, based on self- and clinician-report. 
Five patients were additionally taking an antidepressant or mood- 
stabilizing medication. There were no known cases of extended un- 
treated psychosis. Symptom severity was assessed using the Positive 
and Negative Syndrome Scale (PANSS, Kay et al., 1987), the Beck De- 
pression Inventory-II (Beck et al., 1961 ), and the Beck Anxiety Inven- 
tory (Beck and Steer, 1990) (Table 1 ). Based on rating scores, patients 
could generally be classified as mildly to moderately ill (Leucht et al., 
2005). Patients also endorsed a lower overall quality of life, as as- 
sessed by the abbreviated Quality of Life Enjoyment and Satisfaction 
Questionnaire (Q-LES-Q-l 8, Endicottet al., 1993; Ritsner etal., 2005). 

This study was approved and undertaken in accordance with regu- 
lations established by the Mental Health Research and Ethics Commit- 
tee of Melbourne Health; all participants provided written informed 
consent prior to participation. 

2.2. Magnetic resonance imaging 

Data acquisition was conducted using a Siemens Trio 3 T mag- 
netic resonance scanner at the Murdoch Childrens Research Institute 
(Melbourne, Australia). 

2.2. 1. Resting-state functional MR! 

Over a 1 0-minute period, 250 whole-brain functional MRl volumes 
(36 contiguous axial slices; 3 mm thickness) were acquired while 
participants were awake but performing no explicit cognitive task 
(i.e., standard eyes open resting-state fMRI protocol). Gradient-echo 
echo-planar-imaging parameters were as follows: time to repetition 
(TR) = 2400 ms; time to echo (TE) = 30 ms; flip angle = 90°; field 
of view (FOV) = 210 x 210 mm; in-plane matrix = 64 x 64; and 
in-plane resolution = 3.3 x 3.3 mm. 

2.2.2. Diffusion-weighted MRl 

Diffusion MRl was used to quantify the diffusion of water 
molecules along multiple directions to infer white matter tissue ar- 
chitecture. Sixty diffusion-weighted volumes, each quantifying wa- 
ter diffusion in a unique direction, were acquired using a spin-echo 
echo-planar imaging sequence with the following parameters: b- 
value = 2000 s/mm; 64 consecutive axial slices, 2.0 mm thickness; 
in-plane matrix = 128 x 128; in-plane resolution = 2.0 x 2.0 mm; 
FOV = 256 X 256 mm; TR = 9350 ms; TE = 104 ms; and flip an- 
gle = 90°. Five non-diffusion- weighted volumes were also acquired. 

2.3. Data analysis 

2.3. 1. Preprocessing: resting-state functional MRl 

Preprocessing was performed in MATLAB (version 2012a) using 
the Data Processing Assistant for Resting-State fMRI toolbox (version 
2.2, Chao-Gan and Yu-Feng, 2010). The first 4 volumes (10 s) of the 
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Table 1 

Participant demograpliics. 

SCZ(n = 17) C0N(n = 17) Statistic 



Demograpliics, mean (SD) 



Age, years 


29.9 (6.7) 


27.0 (6.1) 


t32 = 1.31, p = 0.20 


Gender, M/F 


16/1 


1^/4 


— 2.11, p = 0.15 


Education, years 


13.6 (3.0) 


14.6 (1.8) 


— 1 .20, p = 0.24 


Full-scale IQ. 


104(14.9) 


113 (6 1) 


t-,€, — J d'i n — n n77^ 


Handedness, R/L 


15/2 


17/0 


X^ = 2.13, p = 0.15 


Illness and symptom measures, mean (SD) or median [range] 






Age of initial diagnosis, years 


22.5 (4.5) 






Treatment duration, years 


7.5 (4.9) 






PANSS positive 


13.4(6.4) 






PANSS negative 


12.3 (4.5) 






PANSS general 


27.6(6.2) 






Beck Depression Inventory 


16 [0-41] 


2 [0-14] 


Z = 3.61, p< 0.001 


Beck Anxiety Inventory 


7 [0-36] 


2 [0-5] 


Z = 2.83, p = 0.005 


Quality of life'' 


58.4(12.2) 


70.3 (6.7) 


(32 = 3.51, p = 0.001 


Substance use measures, median [range] 








Alcohol, current drinks/month 


2 [0-80] 


9 [0-72] 


Z = 1.86, p = 0.063 


Tobacco, current cigs/day 


0[0-30] 


0[0-15] 


Z = 2.46, p = 0.049 


Cannabis, lifetime uses 


77 [0-3650] 


10 10-700] 


Z= 1.49, p = 0.15 


Ecstasy, lifetime uses 


1 [0-156] 


0[0-30] 


Z=1.50,p = 0.18 


Amphetamines, lifetime uses 


0 [0-365] 


0 [0-4] 


Z = 2.00, p = 0.092 


Cocaine, lifetime uses 


0[0-6] 


0[0-15] 


Z = 0.75, p = 0.63 


Hallucinogens, lifetime uses 


0[0-8] 


0[0-5] 


Z = 0.41, p = 0.76 



Bold entries denote statistical significance; PANSS = Positive and Negative Syndrome Scale. 
"IQnot measured for 1 control and 3 schizophrenia participants. 
'Quality of Life Enjoyment and Satisfaction Scale (abbreviated version). 



acquisition were discarded to allow stabilization of scanner magneti- 
zation and subject orientation. Tiie 36 slices within each whole-brain 
volume were temporally aligned prior to spatial rigid-body coregis- 
tration of the 250 volumes. The realigned data were then normalized 
to IVIontreal Neurological Institute (IVINI) standardized space by first 
estimating a 12-parameter affine transformation mapping the data 
volumes to a standard template image. Nonlinear deformations were 
then estimated to minimize the residual squared difference between 
the data and the template. The warped data were interpolated to a 
spatial resolution of 3 mm isotropic, spatially smoothed using a Gaus- 
sian kernel having a full-width at half-maximum (FWHIVl) of 6 mm, 
and bandpass filtered to isolate relevant low-frequency signal vari- 
ance in each voxel (0.01-0.08 Hz). The following nuisance parameters 
were then regressed out of the data at each brain voxel: six standard 
head motion parameters (3 planes of translation, 3 axes of rotation), 
frame displacement greater than 0.5 mm (including the preceding and 
the two subsequent volumes; note that deleting ( i.e., scrubbing,Power 
et al., 2012) the affected frames yielded very similar results), global 
signal (Keller et al., 2013; Power et al., 2014), and signals from the 
white matter and cerebrospinal fluid. Further exploration indicated 
that the average number of regressed volumes was small and similar 
across groups [tiy gi = 1.98, p > 0.05; average 0.02% (range 0-0.1%) 
in controls; 0.09% (0-0.5%) in patients]. 

The automated anatomical labelling (AAL, Tzourio-lVlazoyer et al., 
2002 ) template was employed to divide the brain into 90 standardized 
regions-of-interest. The cerebellum was excluded from analysis due 
to inconsistent sampling across participants during data acquisition. 
The residuals from the above-described regression were spatially av- 
eraged over the voxels comprising each region. This yielded a set of 
90 regional time series representing 1 0 min of intrinsic (resting-state) 
neural activity for each region. The functional connectivity between 
each pair of regions was quantified independently with the Pearson 
correlation coefficient, resulting in a symmetric 90 by 90 connectivity 
matrix for each participant. 



2.3.2. Preprocessing: diffusion weighted MRI 

Preprocessing was performed using FDT (v2.0), the diffusion 
toolkit in FSL (FIVlRlB's Software Library; http://www.fmrib.ox.ac.uk/ 



fsl). Additional preprocessing steps recommended in the recent lit- 
erature were also performed, as described below (Jones and Cercig- 
nani, 2010). Diffusion-weighted volumes were manually inspected 
for large artefacts and then corrected for stretches, shears, and head 
motion by an affine registration to the first non-diffusion-weighted 
volume using the eddycorrect algorithm. A fractional anisotropy (FA) 
volume was estimated using weighted linear least squares fitted to 
the log-transformed diffusion data. The FA volume was then regis- 
tered to IVlNl space using the FSL nonlinear registration tool FNIRT. 
The nonlinear warp was initialized with an affine registration gener- 
ated with the FSL linear registration tool FLIRT (Jenkinson and Smith, 
2001 ). The registration target was the average of 58 FA volumes from 
healthy adults that were well-aligned to MN1152 space. The inverse 
of the estimated warp field was then applied to the AAL template, 
producing binary masks for the 90 cortical and subcortical regions in 
native diffusion space. 

2.3.3. Probabitistic tractography 

Probabilistic streamline tractography was performed using the 
FSL algorithm ProbtrackX, build 419 (Euler integration; step length: 
0.5 mm; curvature threshold: 80°). Local fibre orientations were 
first estimated using a validated crossing-fibre model (Behrens et 
al., 2003). To map the connectivity between a pair of regions, trac- 
tography was performed with one of the regions serving as the seed 
mask and the other region serving as the termination mask. A to- 
tal of 500 streamlines were initiated from random locations within 
each voxel comprising the seed region. Tractography was then re- 
peated with the seed and termination masks interchanged. The result 
of this was two volumes (connectivity distributions) in which each 
voxel stored the total number of streamlines that it was intersected 
by. These two volumes were averaged voxel-wise and thresholded to 
yield a binary mask of the volume encapsulated by the putative white 
matter pathway interconnecting the pair of regions. The thresholding 
was performed such that any voxel intersected by 1 0 or more stream- 
lines in the averaged map was set to one, while all other voxels were 
set to zero. FA values of all voxels comprising this binary mask were 
averaged to yield a connection-specific measure of white matter in- 
tegrity. This tractography process and tract-specific averaging of FA 
was repeated identically for each participant in all pairs of regions 
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encompassing the functional networks eliciting a between-group dif- 
ference. The FA was adopted as an estimate of structural connectivity 
because of its established link with neurobiological measures of white 
matter integrity (Pierpaoli and Basser, 1996; Kochunov et al., 2007). 
Furthermore, relative to streamlines counts, FA is more normally dis- 
tributed and resides on a continuous scale, thus better justifying the 
use of standard statistical approaches. 

2.3.4. Statistical analysis: functional connectivity 

The network-based statistic (Zalesky et al., 2010) was used to 
identify networks comprising pairs of regions for which the extent 
of functional connectivity significantly differed between the patient 
and control groups. In brief, a two-sample t-statistic was calculated for 
each pair of regions to test the null hypothesis of equality between 
the patient and control groups in mean functional connectivity. All 
possible, (90 x 89) / 2 = 4005, pairs of regions were tested. Among 
the connections exceeding a test-statistic threshold of 3, a search was 
performed to identify any connected network components. A thresh- 
old of 3 was chosen on the basis that it was found to yield consistent 
findings across the parcellation atlases considered and approximately 
corresponds to the nominal uncorrected p-value of 0.005 (see Zalesky 
et al., 2012 for discussion about the choice of threshold). The size of 
each identified network component was measured by the number of 
supra-threshold connections it comprised. Permutation testing was 
then performed to estimate a corrected p-value for each network. 
For each of 10,000 permutations, the size of the largest network was 
recorded to generate a null distribution for calculation of family-wise 
error corrected statistical thresholds. A corrected p-value for a net- 
work identified in the actual data was then estimated by the propor- 
tion of permutations for which a network of equal or greater size was 
identified. 

The AAL template is one of the most prevalently used regional 
parcellations for mapping both the functional and the structural hu- 
man connectome (Zalesky et al., 2010). We utilized the AAL to max- 
imize the number of existing studies with which our results could 
be directly compared without the need to determine a template- 
to-template mapping between discordant regional definitions. Nev- 
ertheless, to evaluate the sensitivity of our findings to the regional 
definitions associated with the AAL, we replicated the above analyses 
using an alternative regional parcellation scheme comprising 1 96 cor- 
tical and subcortical regions of approximately equal volume (Zalesky 
et al., 2010). 

Results were visualized using the software BrainNet Viewer (Xia 
et al., 2013). 

2.3.5. Statistical analysis: structural connectivity and structure- 
function relationship 

Three networks comprising pairs of regions for which the extent of 
functional connectivity significantly differed between the patient and 
control groups were identified (see Results section). For each of these 
disrupted networks, a network-wide measure of white matter in- 
tegrity was calculated. This first involved estimating an edge-specific 
FA value for each constituent connection (based on tractography- 
delineated masks; see Preprocessing: diffusion weighted data) and 
then averaging these connection-specific FA values across all connec- 
tions comprising the disrupted network to yield a single network- 
wide index quantifying white matter integrity (Fig. 1 for examples 
of tractography-delineated white matter pathways). For each of the 
three networks, a two-sample one-tailed t-test was used to test for 
a reduction in network-wide white matter integrity in the patient 
group. This enabled assessment of the hypothesis of co-localization 
of functional and structural connectivity deficits in the schizophrenia 
connectome. 

For each disrupted network, we further tested for a linear rela- 
tionship among participants between their levels of functional and 



Controls Schizophrenia 



a. Right superior frontal cortex - right pallidum 




b. Right superior frontal cortex - left pallidum 




c. Right orbitofrontal cortex - left medial cingulate cortex 




e. Right superior frontal cortex - right thalamus 




Fig. 1. Within each functional network of interest, the average fractional anisotropy 
(FA) between constituent region-pairs was estimated to provide connection and 
networl<-level indices of structural connectivity for each subject. Depicted here are 
examples of white-matter tractography between 5 region-pairs (relevant to Network 
1 ) in healthy controls and patients with schizophrenia. Visual inspection of tractogra- 
phy outputs was important to exclude problems in fibre tracking. However, our analysis 
focused on network-level differences and do not allow to declare changes at the level 
of single connections. Results from confirmatory analyses at the resolution of single 
connections can be found in Supplementary Fig. 2 (a = #20, b — #19, c = #1, d = #6, 
e = #22). 
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Structural connectivity. Functional connectivity estimates were aver- 
aged over all connections comprising each disrupted network to yield 
network-wide functional connectivity averages, comparable to de- 
termination of network-wide structural connectivity indices (above). 
Evidence of a linear relationship in these measures of functional and 
structural connectivity was then assessed for the patient and control 
groups separately using the Pearson correlation coefficient. Fisher's 
method (Fisher, 1 92 1 ) was used to test for between-group differences 
in the strength of these linear structure-function relationships in each 
network of interest. 



3. Results 

3. 1. Functional connectivity 

Compared to controls, patients with schizophrenia showed al- 
tered functional connectivity in three spatially distinct large-scale 
brain networks (Fig. 2). Reduced functional connectivity was evident 
in patients in an anterior functional brain network (pfwe = 0.011; 
Network 1 in Fig. 2). This network comprised 24 brain regions and 
29 connections encompassing frontal, fronto-temporal, and fronto- 
subcortical interactions (Supplementary Table 1). Although changes 
in functional connectivity involved both brain hemispheres, this net- 
work was weighted towards the right hemisphere (Fig. 2). On average, 
the valence of this network was positive in controls (rmean = 0.11; 
ti6 = 6.1; p < 0.001) and negative in patients (rmean = -0.09; 
ti6 = -6.6; p < 0.001). 

Patients with schizophrenia also showed reduced functional con- 
nectivity in a posteriorly localized network (pfwe = 0.017; Net- 
work 2 in Fig. 2). This network consisted of 24 connections be- 
tween 18 regions comprising bilateral fronto-occipital, occipito- 
temporal, occipito-parietal, and parieto-temporal connections (Fig. 
2 and Supplementary Table 1 ). This difference was explained by a 
lack of significant functional connectivity in patients (rmean = -0.03; 
ti6 = -1.19; p = 0.25) while positive connectivity was detected in 
controls (rmean = 0.21; fi6 = 5.84; p < 0.001). 

Finally, patients showed increased functional connectivity com- 
pared to controls in a network primarily encompassing temporo- 
parietal regions (pfwe = 0.007; Network 3 in Fig. 2). This network 
encompassed 33 connections between 30 regions and was weighted 
towards the left hemisphere (Fig. 2 and Supplementally Table 1 ). Re- 
gions encompassing this network were significantly connected in 
patients (rmean — 

0.18; ti6 = 11.2; p < 0.001) but not in controls 
(rmean = -0.05; fig = -1.78; p = 0.09). 

These findings were replicated using an alternative, finer-grained 
parcellation template comprising 192 cortical and subcortical regions 
(Supplementally Fig. 1). 

It is important to remark that the NBS is a method for performing 
network level inference, and thus the null hypothesis was rejected for 
each of the three networks as a whole (p < 0.05), but not for any spe- 
cific connection comprising any network. In other words, between- 
group differences cannot be declared at the resolution of individual 
connections in isolation. 



3.2. Structural connectivity 

In the first network (Network 1 in Fig. 2), patients showed a trend- 
level reduction in white matter integrity (t32 = 1 .63, p = 0.057, Cohen's 
d = 0.59; Fig. 3 and Supplementary Fig. 2). On the other hand, com- 
parable between-group white matter integrity was found in Network 
2 (t32 = 0.57, p = 0.29, Cohen's d = 0.20). Network 3 displayed the 
same pattern of results as Network 1, namely a trend-level reduction 
in white matter integrity in the patient group (f32 = 1.63, p = 0.057, 
Cohen's d = 0.53). 



3.3. Structure-function relationships 

Compared to patients, controls showed a larger correlation be- 
tween the average functional connectivity and white matter integrity 
in Network 1 (Z= 1 .64, p = 0.044). As depicted in Fig. 3, this difference 
was driven by a positive function-structure relationship in the con- 
trol group (r = 0.43, p = 0.032), alongside an absence of correlation 
in the patient group (r = -0.17). These relationships were corrob- 
orated using non-parametric Spearman correlations to account for 
the impact of any outliers in the data (control rho = 0.40; patient 
rho = -0.13). No between-group differences in function-structure 
relationships were found for the second (Z = 0.06, p = 0.48) or the 
third networks (Z= -1.48,p = 0.93). Confirmatory analyses excluded 
that these results were driven by major changes in isolated pairwise 
connections (Supplementary Fig. 2). 

Additional analyses were also performed to account for the po- 
tential confounding influence of age, gender, IQ, or substance use 
(detailed inTable 1 ) on between-group differences. Each variable was 
included as an additional linear regressor in the NBS model assessing 
functional connectivity alterations and in the linear regression assess- 
ing white-matter differences; partial correlations were performed to 
account for these effects in structure-function relationships. Each 
confound was assessed in an independent model. Results indicated 
that the reported connectivity differences were not explained by vari- 
ance related to any of these factors. Similarly, metrics of sub-clinical 
anxiety and depression were not predictive of connectivity or connec- 
tivity relationships in either cohort. Additionally, measures of func- 
tion, structure, and function-structure connectivity coupling in the 
patient group were not significantly correlated with clinical mea- 
sures, including PANSS scores, disease chronicity, and medication 
status. Changes in connectivity obsereed in patients may therefore 
point to trait-like, enduring disruptions governing brain structure and 
function, as opposed to state-like metrics of transient abnormality. 

4. Discussion 

This study sought to assess the relationship between abnormali- 
ties in functional connectivity in schizophrenia and their structural 
basis. We first identified alterations in functional connectivity across 
the whole brain using a cross-sectional experimental design in a 
well-characterized sample of patients with chronic schizophrenia 
and matched controls. The integrity of the white matter pathways 
interconnecting the regions comprising the aberrant functional net- 
works, and the relationship between the observed changes in func- 
tional connectivity and white matter integrity, were subsequently 
assessed. We found that alterations in functional interactions related 
to schizophrenia can be mapped onto both structural deficits and ab- 
normal structure-function relationships in connections involving the 
fronto-striatal, fronto-temporal, and fronto-thalamic regions. Aber- 
rant patterns of functional connectivity were also identified between 
occipito-parietal, occipito-temporal, and parieto-temporal connec- 
tions, but function-structure relationships were relatively preserved 
in these networks. Our findings confirm the utility of a multi- 
modal structural-functional approach to schizophrenia connectomics 
(Skudlarski et al., 2010; van den Heuvel et al., 2013). 

A number of recent findings from theoretical and empirical studies 
have shown that structural connectivity constrains patterns of func- 
tional connectivity (Honey et al., 2009; Hermundstad et al., 2013; 
Honey et al., 2007; Hagmann et al., 2010). In line with this principle 
of brain organization, we have shown co-localized alterations of func- 
tional and structural connectivity in two large-scale brain networks. 
The first network (Network 1, Fig. 2) encompassed the frontal, tem- 
poral, thalamic, and striatal regions, while the second network (Net- 
work 3, Fig. 2) included regions of the parietal and temporal cortices. 
Importantly, in the first network, reductions in both functional con- 
nectivity and structural connectivity in the patient group occurred in 
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Network 1. Control > Schizophrenia 

( 




Network 2. Control > Schizophrenia 




Network 3. Schizophrenia > Control 





Fig. 2. Patients with scliizoplirenia siiowed decreased functional connectivity in two distinct brain networl<s (Networl< 1 and Networl< 2). Networl< 1 encompassed pairwise 
connectivity between frontal, fronto-temporal, and fronto-striatal regions. Network 2 comprised occipito-parietal, occipito-temporal, and parietal-temporal interactions. In a third 
network (Network 3), encompassing parieto-temporal connections, patients showed increased functional connectivity compared to healthy controls. 



the context of a loss of the normal relationship, relative to the control 
group, between function and structure (Fig. 3 ). Similar deficits in func- 
tional and structural connectivity between frontal (Cole et al., 201 1 ), 
fronto-striatal (Fusar-Poli et al., 2011; Fornito et a!., 2013), fronto- 
thalamic (Woodward et al., 2012; Anticevic et al., 2013; Klingner et 
al., 2014), and fronto-temporal (van den Heuvel et al., 2010; Baker et 
al., 2014; Khadka et al., 2013) regions have been consistently found 
in schizophrenia. Specifically, recent research suggests that altered 
fronto-striatal coupling may be reliable risk phenotypes for psychotic 
disorders (Fusar-Poli et al., 2011; Fornito et al., 2013; Fusar-Poli et al., 
2010). Our results add to previous findings by suggesting that altered 
dynamics among regions encompassing Network 1 may result from 
a reduction of the constraints that white matter connectivity pose on 
brain function. 

In contrast to Network 1, reduced structural integrity in temporo- 
parietal connections encompassing Network 3 was paralleled by en- 
hanced functional connectivity in schizophrenia. This result is con- 
sistent with recent findings suggesting that the impact of minor re- 
ductions in structural connectivity can be mitigated by an increase 
in coupling strength among regions comprising the affected network, 



without significantly changing the relationship between functional 
and structural connectivity (Cabral et al., 2012). However, temporo- 
parietal cortices encompassing Network 3 were not functionally inte- 
grated in controls. The observed increase in functional integration be- 
tween temporal and parietal cortices observed in patients may there- 
fore represent a compensatory mechanism induced by more funda- 
mental alterations in the schizophrenia connectome, such as altered 
integration between regions encompassing Network 1 (Fornito et al., 
2013). 

Patients with schizophrenia were also characterized by a lack 
of functional integration between temporal, parietal, and occipital 
cortices (Network 2). In this network, impaired functional connec- 
tivity was found in isolation of significant alterations in the struc- 
tural integrity of the underlying anatomical pathways and changes 
in function-structure relationships. Thus, decreased functional in- 
tegration between regions encompassing Network 2 in patients 
with schizophrenia might be driven by abnormal synaptic plasticity 
(Stephan et al., 2009), as opposed to aberrant anatomical connectivity. 
Accordingly, our findings suggest the existence of common (Network 
1) and distinct (Network 2) pathoaetiologies underlying functional 
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Fig. 3. A significant breakdown in the normal relationship between functional and 
structural brain connectivity was identified only in Network 1. A significant positive 
association between functional and structural connectivity was found in healthy par- 
ticipants, while this relationship was lost in patients with schizophrenia. In addition, 
in Network 1 , patients with schizophrenia showed a trend-level reduction in fractional 
anisotropy (FA) compared to healthy controls ('p — 0.057; Cohen's d — 0.59). 
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and structural changes in the schizophrenia connectome. 

A distinct strength of this study was the recruitment of a cohort 
of individuals suffering from established schizophrenia in the ab- 
sence of severe symptomatology, significant drug use history, or co- 
morbid psychiatric conditions. As such, many common confounds to 
schizophrenia research have been mitigated, increasing confidence 
that the reported results are related to the phenomenology of the 
disorder. However, the modest number of patients recruited under 
such stringent criteria did eschew statistical power. This issue was 
particularly relevant when relating experimental variables to symp- 
tom profiles, given both the small number of samples and the low 
variability inherent to the clinical measures (Table 1 ). Relating altered 
function-structure connectivity coupling to symptoms of schizophre- 
nia therefore remains a compelling avenue for future investigation. 
Replication studies will be important in establishing reliability and 
determining the extent to which these findings relate to different 
symptom profiles. In particular, larger sample sizes may uncover 
more subtle alterations in the schizophrenia connectome and allow 
inference at more conservative statistical thresholds that account for 
multiple comparisons in structure-function relationships across net- 
works or individual connections. Although we found no relationships 
between changes in function, structure, and function-structure cou- 
pling and antipsychotic medication (chlorpromazine equivalents) or 
disease chronicity, the potential impact of medication and illness pro- 
gression on indices of connectivity and their inter-relationships re- 
mains an outstanding question. Similar deficits of brain connectivity 
to those reported herein have, however, been recently reported in 
larger samples of non-medicated individuals who have a genetic li- 
ability for psychosis, as well as in first-episode psychosis (Fornito et 
al., 2013; Crossley et al., 2009). 

In summary, this study highlights the importance of considering 
not only the separable impact of functional and structural connectiv- 
ity deficits in relation to the pathophysiology of schizophrenia, but 
also the implications of the complex nature of their interaction. We 
report a convergence of reduced functional connectivity, impaired 
white matter integrity, and a loss of normal function-structure rela- 
tionships in a network encompassing frontal, striatal, thalamic, and 
temporal regions. Together, these findings suggest that anatomical 
connectivity has less of an influence in constraining functional inter- 
actions between these regions in schizophrenia. 
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